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Einstein-aether theory is general relativity coupled to a dynamical unit timelikc 
vector field. A brief review of current theoretical understanding and observational 
constraints on the four coupling parameters of the theory is given. 
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1. Introduction 

m 

In general relativity (GR), spacetime structure is determined by a dynam- 
ical metric tensor field g a b and nothing else, and the theory is both diffeo- 
morphism invariant and locally Lorentz invariant. Einstein- aether theory 
is the extension of GR that incorporates a dynamical unit timelike vector 
field u a — the "aether" — which breaks the local Lorentz symmetry down to 
a 3d rotation subgroup. Direct coupling of matter to the aether would vio- 



Vh ' late local Lorentz symmetry yet preserve diffeomorphism invariance. This 

paper presents a brief overview of the current theoretical and observational 
status of this theory, assuming that matter does not couple directly to the 
aether. 

The action involving metric and aether is highly constrained. Besides the 
cosmological constant term, the only independent diffeomorphism invariant 
local terms containing no more than two derivatives are 

f 



S ~vtc] ^ {R + K ^aU m V b u n ) d 4 x, (1) 

where R is the Ricci scalar, K^ n is defined as 

KL = c l9 ah g mn + c 2 5<;j h n + c 3 5«5l + c 4 u a u b g mn (2) 

with dimcnsionless coupling constants Cj, and the unit timelike constraint 
on the aether is implicit. (The metric signature is (H ) and the speed of 
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light defined by the metric g a b is unity.) Higher derivatives would be sup- 
pressed by powers of a (presumably) small length, e.g. the Planck length. 
It is assumed here that the aether is aligned at large scales with the rest 
frame of the microwave background radiation. 

Einstein- aether theory — "ae-theory" for short — is similar to the vector- 
tensor gravity theories studied by Will and Nordvedt, 1 but with the crucial 
difference that the vector field is constrained to have unit norm. This 
constraint eliminates a wrong-sign kinetic term for the length-stretching 
mode, 2 hence gives the theory a chance to be viable. An equivalent theory 
using the tetrad formalism was first studied by Gasperini, 3 and in the above 
form it was introduced by Jacobson and Mattingly 4 

2. Newtonian and post-Newtonian limits 

In the weak-field, slow-motion limit ae-theory reduces to Newtonian 
gravity, 5 with a value of Newton's constant Gn related to the parameter G 
in the action (1) by 

G 

Gn = i To' ^) 

1 - C14/2 

where Cu = c\ + C4. (Similar notation is used below for other additive 
combinations of the Cj.) For any choice of the Ci, all parameterized post- 
Newtonian (PPN) parameters 6 of ae-theory agree with those of GR 7,8 except 
the preferred frame parameters 0:1,2 which are given by 8 

«i = -M±^4 (4 ) 

2ci - c\ + c\ K ' 

ai (ci + 2c 3 -c 4 )(2ci +3c 2 + c 3 + c 4 ) 

2 Ci23(2 - C14) 

(This particular way of expressing a^ was given in Ref. 9. The small q 
form of 02 was first computed in Ref. 10.) 

Observations currently impose the strong constraints ai < 10~ 4 and 
ci2 £4x 10 -7 . 6 Since ae-theory has four free parameters Cj, we may set 
01,2 exactly zero by imposing the conditions 8 

c 2 = (-2c 2 -c 1 c 3 + c2)/3c 1 (6) 

c 4 = -c§/ci. (7) 

With (6,7) satisfied, all the PPN parameters of ae-theory are equivalent to 
those of GR. (The parameters 01.2 can also be set to zero by imposing 
C13 = C14 = 0, but this case is pathological, as discussed in section 8.) 
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3. Homogeneous isotropic cosmology 

Assuming spatial homogeneity and isotropy, u a necessarily coincides with 
the 4- velocity of the isotropic observers, and the aether stress tensor is 
just a certain combination of the Einstein tensor and the stress tensor of 
a perfect fluid with energy density proportional to the inverse square of 
the scale factor, like the curvature term in the Friedman equation. 11,5 The 
latter contribution plays no important cosmological role since the spatial 
curvature is small, while the former renormalizes the gravitational constant 
appearing in the Friedman equation, yielding 5 

G 

G CO smo = Z —< —5 7-777- (o) 

1 + (Cl3 +3c 2 )/2 

Since G COS mo is not the same as Gn the expansion rate of the universe 
differs from what would have been expected in GR with the same matter 
content. The ratio is constrained by the observed primordial 4 He abundance 
to satisfy |G CO smo/GN — 1| < 1/8. 5 When the PPN parameters 0:1,2 are set 
to zero by (6,7), it turns out that G COS mo = Gn, so this nucleosynthesis 
constraint is automatically satisfied. 8 



4. Linearized wave modes 

When linearized about a flat metric and constant aether, a>thcory posesses 
five massless modes for each wave vector: two spin-2, two spin-1, and one 
spin-0 mode. The squared speeds of these modes relative to the aether rest 
frame are given by 

spin-2 1/(1 - C13) (9) 

spin-1 (ci - \c\ + |c§)/cu(l - C13) (10) 

spin-0 ci 23 (2- cu)/cu(l - c 13 )(2 + c 13 + 3c 2 ) (11) 

The corresponding polarization tensors were found in one gauge in Ref. 12 
and in another gauge in Ref. 9. The energy density of the spin-2 modes is 
always positive, while for the spin-1 modes it has the sign of (2ci — c 2 + 
c§)/(l — C13), and for the spin-0 modes it has the sign of Ci4(2 — C14). 13 ' 9 
(These reduce to the results of Ref. 14 in the decoupling limit where gravity 
is turned off.) 

These squared speeds correspond to (frequency/wavenumber) 2 , so must 
be non-negative to avoid imaginary frequency instabilities. They must 
moreover be greater than unity (super-luminal), to avoid the existence of 
vacuum Cerenkov radiation by matter. 2 (The strongest constraints arise 
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from the existence of ultra high energy cosmic rays.) And the mode en- 
ergy densities should be positive, to avoid dynamical instabilities. With the 
«i,2 = conditions (6,7) imposed, all of these conditions are met for all of 
the modes if and only if c± = C\ ± C3 are restricted by the inequalities 8 

< c+ < 1 (12) 

0< c_ <c+/3(l-c+). (13) 

Interestingly if the mode speeds are instead required to be less than unity 
(sub- luminal) , then the spin-1 and spin-0 energy densities are negative. 
Hence not only the Cerenkov constraint, but also energy positivity (together 
with ai.2 = 0) requires mode speeds greater than unity. 

Note that when (7) holds, we have C14 = 2c + c_/(c + + c— ), which sat- 
isfies < C14 < 2 when the constraints (12,13) hold. Thus in particular 
the condition for attractive gravity mentioned in section 2 need not be 
separately imposed, and c\\ is non-negative. 

5. Primordial perturbations 

Given the same Gn, and assuming the PPN parameters 0:1.2 vanish, the pri- 
mordial power in cosmological spin-0 and spin-1 perturbations is unchanged 
relative to GR, while the power in spin-2 perturbations differs from that in 
GR by the factor (1 - c M /2)(l - c^) 1 / 2 . 14 ' 15 When the constraints (12,13) 
arc satisfied this factor is smaller than unity, hence these spin-2 perturba- 
tions are even more difficult to detect than in GR. As for the late time 
evolution of these perturbations, neutrino stresses in the radiation domi- 
nated epoch source the spin-1 mode, which leads to modified matter and 
CMB spectra. The effect is rather small however, and is degenerate with 
matter-galaxy bias and with neutrino masses. 15 

6. Radiation damping and strong self-field effects 

If the fields are weak everywhere (including inside the radiating bodies), 
and the PPN parameters 0:1,2 vanish, radiation is sourced only by the 
quadrupole. Waves of spins 0, 1 and 2 are radiated, and the net power is 



.•;-' 



given by (GnA/5)Q^, where Qij is the quadrupole moment and A — A[ 



c, 



is a function of the coupling parameters Cj that reduces to unity in the case 
of GR. 9 Agreement with the damping rate of GR (confirmed to ~ 0.1% 



in binary pulsar systems 6 ) can be achieved by imposing the condition 
A[ci] — 1, which is consistent with the constraints (12,13). 
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Compact sources with strong internal fields such as neutron stars or 
black holes can be handled 16 using an "effective source" dynamics specified 
by a worldline action integral 

5 = -mo [dT[l + a(v a u a -l) + a'(v a u a -l) 2 + ...}, (14) 



where v a is the 4-velocity of the body, u a is the local background value 
of the aether, and a and a' are constants characterizing the body, called a 
"sensitivity parameters" or just "sensitivities" . The sensitivites scale as Cj 
for small Cj. 

The effects of nonzero sensitivities on two-body dynamics and radiation 
rates lead to a number of phenomena that are constrained by observations, 
including violations of the strong equivalence principle, modifications of the 
post-Newtonian dynamics, modifications of quadrupole sourced radiation, 
and both monopole and dipole sourced radiation. When 0:1,2 = 0, all of 
these constraints are met provided the sensitivities are less than <~ 0.001, 
which will certainly be the case if Cj < 0.01. 16a To be more precise would 
require knowing the actual dependence of the sensitivities on the c, , which 
has so far only been determined for a and only at leading order (where a 
vanishes when 0:1,2 = 0). (The speed V of the observed binaries with respect 
to the background aether frame can be neglected in formulating these con- 
straints provided V < 10~ 2 , which is easily satisfied for any known proper 
motion relative to the rest frame of the microwave background radiation. 16 ) 

7. Spherically symmetric stars and black holes 

Unlike GR, ac-theory has a spherically symmetric mode, corresponding to 
radial tilting of the aether. For each mass, there is a two parameter fam- 
ily of spherically symmetric static vacuum solutions, rather than a unique 
solution as in GR. 18 Asymptotic flatness reduces this to a one parameter 
family 7,18 The solution outside a static star is the unique solution for a 
given mass in which the aether is aligned with the Killing vector. 1 This 
"static aether" vacuum solution depends on the Ci only through the combi- 
nation C14, and was found analytically (up to inversion of a transcendental 
equation). 18 It is stable to linear perturbations under the same conditions 
as for stability of flat spacetime, with the exception of the case C123 = 0. 19 



a This corrects an error in version 1 of Ref. 16, where a is said to scale as c|. (Also the 
a prefactor C14 in Eqn. (70) should be deleted.) As a result of this correction, the likely 
constraints on Ci are an order of magnitude stronger, as stated here. 1 
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The solution inside a fluid star has been found by numerical integra- 
tion, both for constant density 18 and for realistic neutron star equations of 
state. 20 The maximum masses for neutron stars range from about 6 to 15% 
smaller than in GR when C14 = 1, depending on the equation of state. The 
corresponding surface redshifts can be as much as 10% larger than in GR 
for the same mass. Measurements of high gravitational masses or precise 
surface redshifts thus have the potential to yield strong joint constraints on 
C14 and the equation of state. The radius of the innermost stable circular 
orbit (ISCO) differs from the GR value 6G^M by a small term of relative 
order about 0.03ci4. 

For black holes, the condition of regularity at the spin-0 horizon selects 
a unique solution from the one-parameter family for a given mass. 21 When 
a black hole forms from collapse of matter, the spin-0 horizon develops in 
a nonsingular region of spacetime, where the evolution should be regular. 
This motivated the conjecture that collapse produces a black hole with 
nonsingular spin-0 horizon, which has been confirmed for some particular 
examples in numerical simulations of collapse of a scalar field. 22 

The black holes with nonsingular spin-0 horizons are rather close to 
Schwarzschild outside the horizon for a wide range of couplings; for instance, 
the ISCO radius differs by a factor (1 + 0.043ci + 0.061c 2 ), in the case with 
C3 = C4 = and C2 fixed so that the spin-0 speed is unity 23 (This expansion 
is accurate at least when c\ < 0.5. No solution with regular spin-0 horizon 
exists in this case when c\ > 0.8.) Inside the horizon the solutions differ 
more, but like Schwarzschild they contain a spacelike singularity. Black 
hole solutions with singular spin-0 horizons have been studied in Ref. 24. 
These solutions can differ much more outside the horizon. Quasi-normal 
modes of black holes in ae-theory have been investigated in Refs. 25. 



8. Special values of Ci? 

The first case to be examined in detail 26 ' 4 was C13 = c^ = C4 = 0, i.e. 
the "Maxwell action" together with the unit constraint on the vector. The 
PPN result for 02 (5) is infinite in this case, and the spin-0 mode speed is 
zero. The perturbation series used in the PPN analysis is thus evidently 
not applicable. Independently of that however, other problems with this 
case have been identified, such as the formation of shock discontinuities 4 ' 27 
and a possibly related instability 19 

Assuming now that 0:1,2 = and the constraints (12,13) are satisfied, 
and putting aside the case c\ — C3 = which is not covered by existing 
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PPN analyses, all but one of the cases in which one of the Cj vanishes, 
or in which one of C13, C14, or C123 vanishes, have the property that the 
spin-1 mode speed (10) diverges while the energy of that mode is nonzero. 
It seems very unlikely that such cases are observationally viable, although 
they have not been examined carefully. The exception is the special case 
C3 = C4 = 2ci + 3c2 = 0, with 2/3 < c\ < 1. This large value of c\ is 
probably inconsistent with the strong field constraints from orbital binaries, 
but as mentioned above those are not yet precisely known because the 
sensitivity parameters have not yet been computed for neutron stars, so 
this case is not yet ruled out. 

9. Conclusion 

Einstein-aether theory is an intriguing theoretical laboratory in which grav- 
itational effects of possible Lorcntz violation can be meaningfully studied. 
There is a large (order unity) two-parameter space of Einstein-aether the- 
ories for which (i) the PPN parameters are identical to those of GR, (ii) 
the linear perturbations are stable and carry positive energy, (iii) there is 
no vacuum Cerenkov radiation, (iv) the dynamics of the cosmological scale 
factor and perturbations differ little from GR, (v) non-rotating neutron star 
and black hole solutions are close to those of GR, but might be distinguish- 
able with future observations. Radiation damping from binaries, imposes 
an order 0.001 constraint on one combination of the parameters. Strong 
self-field effects in neutron stars and black holes produce violations of the 
strong equivalence principle and higher order post-Newtonian effects which 
will constrain all the parameters Cj to be less than around 0.01, presum- 
ing that the sensitivity parameters for neutron stars (which have not yet 
been computed with the required precision) turn out to have the expected 
magnitude. 
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